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Quantum constraints of the type Q\ipphys) = can be straightforwardly imple- 

^5 ■ mented in cases where Q is a self-adjoint operator for which zero is an eigenvalue. 

^D I In that case, the physical Hilbert space is obtained by projecting onto the kernel 

of Q, i.e. Hphys = ker Q = ker Q* . It is, however, nontrivial to identify and 

Mh, project onto Hpi^yg when zero is not in the point spectrum but instead is in the 

■y \ continuous spectrum of Q, because then ker Q = 0. 

C^ ■ Here, we observe that the topology of the underlying Hilbert space can be 

^ . harmlessly modified in the direction perpendicular to the constraint surface in 

' such a way that Q becomes non-self-adjoint. This procedure then allows us to 

conveniently obtain Hp^yg as the proper Hilbert subspace Hphys = ker Q* on 

^ ■ which one can project as usual. In the simplest case, the necessary change of 

H ! topology amounts to passing from an L^ Hilbert space to a Sobolev space. 



1 Introduction 

Numerous classical dynamical systems are distinguished by the presence of constraints, 
which, in a phase space formulation, act to restrict the system to the constraint hy- 
persurface, a submanifold of the original phase space with a positive co-dimension. 
In practice, this restriction is accomplished by one or more real constraint functions, 
(paip, q) = 0, a = 1,2, . . . , A, A < oo, which are nondynamical — no time derivatives — 
and serve to constrain the system to the constraint hypersurface. When dealing with 
concrete examples it is useful to further categorize a system of constraints (along with 
their associated Hamiltonian) into classes (first and second), as well as other subdi- 
visions (closed, open, irreducible, reducible, regular, irregular, etc.); these categories 
are well described in the literature and are not reviewed here |1| . For purposes of the 
present paper, it is sufficient to focus on the several constraint functions themselves, 
and we need not be too concerned about any specific subclassification of the set of 
constraints. 

Our principal interest lies in quantization, or more particularly, with the quantum 
theory of constraints. We assume that as quantum operators, the constraints $„, 
a = 1, . . . ,A, are represented by self-adjoint operators determined in some fashion 
from the classical constraint functions by some consistent but unspecified quantization 
procedure. Just as the classical constraints act to restrict the system to a subset of the 
original classical phase space, it is the role of the quantum constraints, in like manner, 
to restrict the system to a subset of the original quantum mechanical phase space. 
Since the "quantum mechanical phase space" is a Hilbert space, such a restriction is 
ideally imposed by the several constraint conditions $„ \'^)phys = 0, a = 1, ■ ■ ■ ,A. As 
a linear equation, it follows that the vectors \il')phys form a linear space, and since the 
constraint operators are self adjoint, the given space is closed, hence a subspace Hphys, 
the presumed physical Hilbert space, in the original Hilbert space, H. For evident 
reasons we focus attention on those cases where dim(if . ) > 0. 

The foregoing scenario does indeed hold for certain families of constraints, namely, 
in cases where the constraint operators $„ each have a set of simultaneous eigenvec- 
tors with eigenvalue zero. When that is the case, we may also consider the single 
constraint operator X^ = S^^^$^, regarded as a self-adjoint operator, and observe 
that X^ \'ip)phys = is completely equivalent to the several equations $q, \ip)phys = 0, 
a = 1, . . . ,A. An example of this kind of constraint situation is given by $q, = Ja, 
a = 1,2, 3, where the operators Ja satisfy the Lie algebra for S0(3) [or SU(2)] and the 
condition E J^ \'ip)phys = corresponds to a restriction to spherically symmetric states. 
We also observe that we can also set Hp^yg = EH, where E = E* = lE^ is a projection 
operator, which in the present case is E = E(EJ^ = 0). For comparison purposes with 
what is to come we also note, equivalently, that 

E = E(EJ2 < ch^) = E(SJ2 = 0) , (1) 

where c < 2 for SO (3) [or c < 3/4 for SU(2)]. The discussion and example of the 



present paragraph refer to the ideal situation regarding the constraint operators. 

More generally, the constraint operators do not fulfill the ideal (Dirac) criteria 
given above. In fact, it frequently happens that the set of constraint operators have 
no nonzero eigenvector with eigenvalue zero. This situation may arise in two funda- 
mentally different ways. One of these ways refers to cases where X^ = S^^]^$^ has a 
discrete spectrum (in the vicinity of zero) which does not include zero. An example of 
this situation is given by the two classical constraints (pi = p = and 02 = Q' = (for 
the same degree of freedom), which then become the quantum constraints $1 = P, and 
$2 = Qi two operators which satisfy the Heisenberg commutation relation [Q, P] = i 
(with ^ = 1). The ideal equations, P \ip)phys = and Q \ip)phys = 0, imply (modulo 
domain issues) that [Q,P] \ip)phys = i \4')phys = 0, namely, that \ip)phys = and thus 
Hphys = 0, which is unacceptable. (This is a typical case of second-class constraints.) 
To avoid this situation, we replace the ideal conditions by the choice Hphys = E-f^, 
where in the present case E = E(P^ + Q"^ < h) = |0)(0|, the projection operator onto 
the harmonic oscillator ground state, thus leading to a one-dimensional Hphys- More 
generally, we accommodate this kind of situation by the criterion that Hphys = Ei7, 
where 

E = E(X2 < 5{n)^) (2) 

and 6{h) is an ^-dependent regularization parameter to be chosen on a case- by-case 
basis. The choice of a quadratic combination of constraints is not written in stone, and 
a discussion of alternative choices is presented elsewhere 0]. 

The second manner in which the idealized situation may fail arises when the opera- 
tor X^ = S;^^]^$^ has its zero in the continuous spectrum. An example of this situation 
is given by the single classical constraint = g = 0, which then becomes the quantum 
constraint ^ = Q. Other examples could arise from two degrees of freedom for which, 
in an obvious notation, e.g., (i) $1 = Qi, $2 = Q2, or (ii) $1 = Qi, $2 = P2, or (iii) 
^ = Ql + Ql — I, etc. In these cases X^ = T,^^^^'^ has its zero in the continuum. As 
a consequence E(X^ = 0) = 0, which we deem to be unacceptable. In place of this 
idealized condition we once again choose 

E = E(X2 < 6{n)'') (3) 

where 6{h) > 0. Observe, in the present case, that for any 5 > 0, it follows that 
Hphys = Eif is an infinite- dimensional (regularized) physical Hilbert space. Although 
it is actually possible to work with this regularized space, say in cases where 6 is 
extremely small, e.g., 6 = 10~^°°°, in some natural units, it is analytically preferable 
if we are able to take the limit 6 —>■ 0. However, this limit cannot be taken in any 
straightforward fashion since 

lim(A|E(X2< 52)10)^0 (4) 

0— >u 
for any pair |A), |0) G if , whenever the zero of X^ lies in the continuum. 



The problem of quantum constraints with their zero in the continuum is well known 
and has often been studied in the literature. Briefly summarized, such problems have 
been studied, e.g., by (i) introducing Gel'fand triplets [0], (ii) specialized algebraic 
representations 0], or (iii) rescaled limits within suitable subspaces 0. 

As one version of type (iii) above, we imagine working in a representation in which 
X is diagonalized (as x), and therefore 

(A|E(X2 < 6')\4>) = j f_^ W^(P{x,y) dxda{y) , (5) 

where the variable y corresponds to any degeneracy. In this form it is clear, as 6 —>■ 0, 
that the right-hand side vanishes. To overcome this situation, let us first restrict 
attention to a subset of functions, e.g., 

"^ = {xix,y) : xix,y) G {polynomial (x,?/) e'^'' +^ ^}} , (6) 

and consider X = \q E V and (f) = (po E V , so that 

(Ao|E(X2 < 6')\(Po) = j j\ 'M^Mx,y)dxda{y) . (7) 

In the present case it follows that 



lim(25)-MAo|E(X2<<52)|0o) 



lim(25) W / Xo{x,y)(t)o{x,y)dxd(r{y) 



Ao(O,l/)0o(O,l/)c?a-(y) 
= (Ao, 0o) ■ 

This final expression defines a sesqui- linear form characterizing a pre-Hilbert space. 
Completion of the pre-Hilbert space in the usual fashion [i.e., inclusion of limits of 



Cauchy sequences in the inner-product-induced norm ||0o|| = \/{4>q,4>o)j plus the iden- 
tification of elements as equivalence classes of functions as necessary] leads to the true 
physical Hilbert space, Hphyg, in which the constraint condition X"^ = is finally 
satisfied. 

In the previous discussion it is noteworthy that we chose a representation for the 
constraint operators in order to achieve a successful rescaling and limit as the regular- 
ization was removed. Other schemes may avoid the regularization and its subsequent 
removal, but they all require a representation to be introduced. 

A natural question then arises whether or not it is possible to devise a procedure to 
enforce a quantum constraint whose zero lies in the continuum in an abstract fashion, 
i.e., is it possible to impose such a constraint without introducing any specific represen- 
tation whatsoever. It is this question that we address in the remainder of this paper, 
and it is noteworthy that the answer to that question is in the affirmative. 



The key to obtaining this affirmative answer is worth noting. The traditional view 
of obtaining a physical Hilbert subspace in the case of zero in the continuous spectrum 
is to change the Hilbert space, essentially, by imposing the constraint as a distribution. 
The new method introduced in this paper changes the underlying topology of the ex- 
isting (pre-)Hilbert space and then is able to impose the constraint with a zero in the 
continuum in the same way one imposes a constraint with a discrete zero! 

In this section, we have used the symbols P and Q in their traditional Heisenberg 
and Schrodinger senses. Note well that the same symbols P and Q appear in the 
subsequent sections, but they are used in a much wider sense. In particular, Q is used 
as a generic constraint operator with its zero in the continuum, while P is used as a 
maximally symmetric operator satisfying [Q, P] = il on a suitable domain which need 
not even be dense. 

2 The continuous spectrum and topology 

Let us begin with the observation that, for topological reasons, in the case of zero being 
in the continuous spectrum of Q, there cannot exist a projection from the original 
Hilbert space H to a physical Hilbert subspace Hphys'- 

To see this, let us consider the simplest case, namely where the spectrum of the 
constraint operator Q is nondegenerate, and is given by an interval, say I = [0,cxd). 
In the spectral representation of Q, the original Hilbert space H is therefore the space 
H = L'^{I) of square integrable functions f E H over the interval /. 

We would like to find a physical Hilbert subspace Hphys corresponding to "Q = 0" . 
It should be a one-dimensional vector space of "function values at 0" . We could then 
identify this space with C. To this end, let us consider the linear functional defined 
by the property that it maps every continuous function f E H onto its value at zero, 
i.e.: 

: / ^ /(O) (8) 

Intuitively, one might assume that is the desired projection onto a one-dimensional 
"physical" subspace of function values at zero. 

However, is not a projection. The problem is that, even though we restricted 
the domain of to continuous representatives of square integrable functions, is not a 
continuous map: Consider, for example, the sequence {(7„(x)}^g]jy[- of continuous square 
integrable functions 

g^ix) := A e-^'-\ (9) 

for an arbitrary constant A. In the topology of the Hilbert space H, namely the norm 
topology induced by the L^ scalar product of H , the sequence clearly converges to the 
null vector in H, i.e to the function g[x) = . Namely, lim„_>oofi'n = fi* in the sense 
that hm„^oo \\g — 5'nlP = liiii„^oo X^ c^a^ A*Aexp{—2nx^) = 0. On the other hand, 
lim^^oo 0(fl'n) = A, which proves the discontinuity. 



This means the following: We can take any continuous square integrable function 
f E H and add a null sequence such gn{x)- Then, in the limit we recover the vector 
f E H. However, crucially, as far as the original topology of the Hilbert space H is 
concerned, we can in this way give </>(/) any arbitrary value that we wish, i.e. (f) is 
ill-defined. 

Intuitively, the underlying reason for the ill-definedness of is of course that even 
continuous functions f E H can be arbitrarily sharply peaked over a point, say zero, 
meaning that functions f (z H could possess arbitrarily different values at zero - while 
their L^-distance is arbitrarily small, i.e. while such functions can still be arbitrarily 
close within the topology of H. 

Thus, the situation is that is noncontinuous and therefore also unbounded. Since 
projectors are of course bounded, is not a projection! Another way to see this is that 
by Riesz' theorem, 0, being noncontinuous, is not contained in (the dual of) H, and 
can therefore not be projected onto. 

The case of real physical interest is of course where the spectrum of Q is degenerate 
at zero, and where one therefore expects the physical Hilbert space to be higher dimen- 
sional. But it is clear that the fundamental problem, namely that is not a projector, 
persists to all nontrivial cases - as long as we stick to the original topology of H . Our 
aim is therefore to suitably modify the topology of the original Hilbert space in order 
to be able to project onto the physical subspace. 

3 Hints from the theory of distributions and Sobolev 
spaces 

Before we describe our solution to the general problem, let us recall some basic facts 
from the theory of distributions and Sobolev spaces, facts which can be directly used 
to solve the problem for the simple nondegenerate example which we just discussed. 
Our general method will be motivated by this example. 

We begin by recalling how, in spite of our arguments above, the map can be 
made into a continuous and bounded map, in fact, into a projection. The price to be 
paid is that the topology of the original Hilbert space must be changed. Intuitively, the 
topology of the function space needs to be changed in such a manner that functions 
which differ only by sharp peaks (which make arbitrarily little difference in the L^ 
topology) are separated in the new topology. 

The basic observation underlying the theory of Sobolev spaces is that a sufficiently 
improved topology is induced by a simple new scalar product on the function space, 
see e.g. 0. The new scalar product is arranged to be sensitive to the rate of change of 
functions. In the norm topology induced by this scalar product even those functions 
can be separated which differ only by spikes that are so sharp that the L^ topology 
cannot distinguish them (e.g. those differing by sequences such as the sequence {gn{.x)} 



ofEq.|). 

In particular, while the scalar product of H is of course 



(/l|/2)= dxn{x)h{x) (10) 



/ 



the scalar product of the Sobolev Hilbert space H^ is defined, using distributional 
derivatives, as: 

/r(x)/,(x) + ^^*^^)^^^(^)" 



Ui\f2) = j^dx 



dx dx 



(11) 



Note that we will use throughout the notation (|) for the scalar product in the initial 
Hilbert space and (|) for the new scalar product. 

Let us check whether the functional is a continuous functional over H^\ To 
this end, we recall, by Riesz' representation theorem, that continuous functionals over 
a Hilbert space can be identified with the Hilbert space vectors themselves, via the 
scalar product action. Thus, if (f) is indeed a continuous functional over if ^, then we 
should now be able to identify a representation of as a vector in HH In fact, 

0(x) = e"^ (12) 

is the representation of in H^. It is normalized with respect to the scalar product 



Eq.ll of H^. The reader may check in a short calculation that indeed (</>,/) = /(O) 
for all / G H^. Thus, 0, which maps continuous functions onto their values at a point, 
or more precisely 

H := 10) ® (01 (13) 

becomes a projection, as desired. It is also clear that this could not have been achieved 
without a change in topology. 

Remark: Higher Sobolev spaces 

We should also mention that our particular choice of change of topology - to 
move from H to H^ - is not the only possible one to achieve our goal of making a 
projection: To see this, let us be more precise about the properties of the functions in 
H^ as opposed to those in H . For a function f{x) to be in H^ it must not only be 
square integrable: f{x) must also be the indefinite integral of some square integrable 
function ^^df{x)/dx". We recall that a function f{x) is an indefinite integral exactly if 
it is absolutely continuous (on any countable set of non-overlapping intervals). Thus, 
moving from the Hilbert space H = L^ to the Hilbert space H^ improves the continuity 
of the functions f{x) to absolute continuity. However, the derivatives df{x)/dx of 
functions in H^ need only be in L^, not H^. Thus, in H^ it is still not possible to 
project df{x)/dx onto its value at a fixed point. In fact, df{x)/dx need only be defined 
almost everywhere. As we are here only interested in projecting onto the values of 
functions and not their derivatives this is not a problem for our purposes. 



Nevertheless, for completeness, let us recall how the topology could be arranged 
to improve also the behavior of the functions' derivatives. As the theory of Sobolev 
spaces shows, it is sufficient to this end to use higher derivative operators in the scalar 

product, i.e. to use e.g. the scalar product (/i, /2) := Jjdx Yl]=o 2cj jilf for some 
n > 1 and its induced norm topology. In this way one obtains the higher Sobolev 
spaces if". Obviously, H D H^ D H^ D H^... . In the function space ii", we can 
project the derivatives d"^~^ f / dx"^~^ of functions f{x), up to m = n — 1, onto their 
values at any given point. It is clear that, while we could use any higher Sobolev space, 
H^ does suffice for our purposes here. 

4 The new method 

Let us now consider the general problem of projecting onto the physical subspace in 
the case of zero being in the continuous spectrum of a self-adjoint constraint operator 
Q whose spectrum may be arbitrarily degenerate. We would like to generalize the 
procedure used when zero is only in the point and not in the continuous spectrum of 
the self-adjoint constraint operator Q, in which case one can straightforwardly define 
the physical subspace as the kernel of Q, i.e. 

Hphys ■■= ker Q = ker Q* . (14) 

In the case where zero is only in the continuous and not in the point spectrum of Q, the 
problem is of course that then keiQ = kevQ* = 0, which is obviously not the desired 
subspace. 

Our main idea in this paper is therefore to treat this case by modifying the topology 
of the initial Hilbert space H in the direction away from the "constraint surface" in a 
manner analogous to passing from an L^ Hilbert space to the Sobolev space H^. (In 
the general case, where zero is also in the point spectrum, one may project out the 
kernel of Q as usual, before applying our procedure). 

The new Hilbert space, which we will call H, will be a subspace of the original 
Hilbert space when considered as a vector space, and the action of all operators, their 
commutation relations, etc, therefore remains unchanged. Also Q is of course still the 
same, as a linear map. Crucially, however, H will be different as a Hilbert space, being 
equipped with a new scalar product. This scalar product changes the induced norm 
topology in the direction of the constraint degree of freedom. As a consequence, the * 
structure will change and the operator Q will no longer be self-adjoint, i.e. Q ^ Q*\ 
As our main finding, we will show that this modification enables us to identify the 
physical subspace as: 

Hphys ■= ker Q* (15) 

Since Hphys is therefore a proper Hilbert subspace, namely Hphys C H, this means that 
once one has passed to the new topology, i.e. from H to H, one can again implement 

8 



the constraint by projection. 

Explicitly, in order to carry out the program of suitably modifying the topology in the 
direction away from the constraint surface, we begin by completing the description of 
the degree of freedom which is to be constrained: Namely, in addition to the constraint 
operator Q, let us also consider a variable P which is conjugate to Q, i.e. a maximal 
symmetric operator P which obeys: 

[Q,P]=i\ (16) 

By definition, therefore, the domain of P is the maximal domain Dp & H on which 
this commutation relation holds and on which P is symmetric. We remark that while 
in general P will not be self-adjoint on Dp - in particular it never is if Q is positive - 
the property of being maximally symmetric will suffice. As we will see, it will also not 
matter for our purposes that the domain Dp G H on which the commutation relation 
Eq.|TB| holds will in general not even be dense in the original Hilbert space if - as it 
clearly will not be if Q also possesses a point spectrum. Intuitively, this is because for 
our purposes only the part of the spectrum around zero matters. 

In analogy with the Sobolev space example for functions, let us now consider the 
domain Dp* and let us change the topology on it to obtain a new Hilbert space H: 
Namely, as a vector space, we define H to be identical to Dp* while we equip Dp* with 
a new scalar product which then induces a new norm topology. Denoting the scalar 
product in the original Hilbert space H, and thus on Dp*, by ( | ), we define the scalar 
product { \ ) on H through: 

{V1\V2) := {Vi\v2) + {P*Vl\P*V2) (17) 

In other words, H is the graph Hilbert space of P*. In the simple case where the 
spectrum of Q is purely continuous and nondegenerate this Hilbert space is the Sobolev 
space H = H^. 

An advantage of our functional analytic definition is that, unlike in our discussion 
of the simple example of the nondegenerate spectrum and its Sobolev space, we now 
no longer need to work in the spectral representation of Q. Let us denote the domain 
of Q in if by Dq (consisting of all vectors (j) E H for which Q(f) has finite norm with 
respect to the new scalar product given in Eq . p^) . On Dq the operator Q is no longer 
self-adjoint, Q y^ Q* (and not even symmetric). We then define the physical Hilbert 
subspace Hphyg C ii as the kernel of Q*: 

Hphys := ker Q* (18) 

This definition has a simple interpretation: We require that if we act with Q on any 
vector in its domain in H then the resulting vector is orthogonal to all physical vectors. 
In other words, we define the physical subspace Hphys as the orthogonal complement 
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of the range of Q in the Hilbert space H. 

For added clarity, let us be fully precise regarding the definition of Hphys through 
Eq.|l8|. As a vector space, Hphys is given by 

Hphys ■={\v)eH = Dp,\\/ we Dq : {Qw\v) + {P*Qw\P*v) = o} (19) 

i.e. all its vectors are also vectors in H, when considered as a vector space. The scalar 
product in Hphys is given by ( | ) in Eq.p!7|. 

To summarize, we begin by completing the picture of the degree of freedom which 
is to be constrained, namely by augmenting Q by a symmetric operator P obeying the 
commutation relation [Q, P] = il on its maximal domain Dp in the original Hilbert 
space H. Second, we change the scalar product and consequently the induced topology 
on the domain Dp* G H to obtain the graph Hilbert space, H, of P* with the scalar 
product Eq.|Tl|. Third, we identify the physical subspace Hphys as the proper Hilbert 
subspace Hphys = ker Q* of the Hilbert space H. Thus, after passing from H to H we 
can again project onto a physical subspace which is a proper Hilbert subspace. 

Remark: The special case of positive Q 

If Q is a positive self-adjoint operator with zero being in the purely continuous 
spectrum, i.e. when zero is actually a boundary of a piece of the continuous spectrum, 
then there is additional structure which can be used: 

Namely, consider in this case Dp* modulo Dp in the new topology, i.e. consider 
Hb '■= Dp* Dp. By von Neumann's theory of self-adjoint extensions of symmetric 
operators, the domains of P* and P differ exactly by the space of boundary functionals, 
and we have the von Neumann formula 

Dp* Q Dp = N^ ® N- (20) 

where N^ are the deficiency spaces of P. Since we assume that zero is a boundary of 
the spectrum we can therefore conclude that Hphys C A^"*" © A^~. As a consequence, we 
can say that all physical vectors must be a linear combination of vectors in the kernels 
of either the operator (P* + i) or (P* — i). In concrete representations, this fact can 
yield useful differential equations. 

5 Examples 

Let us illustrate the working of the new method with simple examples in non relativistic 
quantum mechanics. We will choose constraint operators Q which are functions only 
of the position operators Q = (^(x). In this way, the associated classical constraint 
surface, namely the zero set of Q = read as a classical equation, will be obvious. 
We will show that our new method projects on the Hilbert space of functions over the 
constraint manifold. 

10 



5.1 Example: Projecting onto a point on the line 

As a first example, let us reconsider the simple case of a constraint Q = for a self- 
adjoint operator Q whose spectrum is continuous at zero and nondegenerate. The 
nondegeneracy means of course that we expect the physical Hilbert space to be one- 
dimensional. This example will basically be the same as the one in which we discussed 
Sobolev spaces. In order make this case slightly nontrivial let us choose the spectrum 
of Q as: 

/=[0,1] U {2,3,...} (21) 

Thus, in the spectral representation of Q, our starting Hilbert space H is H = 
L^(0, 1) © P, with scalar product: 

(0i|02) = / dq<Pl{q)Mq) + Y.<P»Mn) (22) 

We begin by considering a maximally symmetric operator which obeys [Q, -P] = «1, 
namely P = —idq on its domain Dp. Functions in Dp are absolutely continuous, 
square integrable, vanish at the interval boundary, and their almost everywhere defined 
derivative is also square integrable. We note that, while Dp is not dense in H - since it 
lacks the eigenspaces to the discrete eigenvalues - this will not matter when we project 
onto "g = 0". 

We now equip Dp with the new scalar product, Eq.^, to obtain H: 

(0l|02) := (0l|02) + (P*0l|P*02) (23) 

dq {<Pliq)Mq) + id,Mq)rdMq)} (24) 



Next, we define the (Sobolev) Hilbert space H as the graph Hilbert space of P*, i.e. 
as the domain of P*, with the scalar product Eg .[2^. 

Now we are ready to calculate the physical subspace Hp^ys = ker Q*. Explicitly, 
according to Eq.|l^, the condition for vectors ip ^ H to he in the physical domain Hp^ys 
now reads 

[' dq {<P*{q)q^{q) + (d^* (q)) d^q)} =0 ^ (f) E Dq (25) 

JO 

This yields the condition, again for all (p G Dq, 

I'dq {r{q)qHq) " q<P*{q)d',i^{q)} + W iq)dMl)]l = (26) 

which means that 

g(5j - l)^(g) = (27) 

and 

5,V^(g)l,=i = 0. (28) 

11 



There is only one solution to Eqs.^^ up to normalization: 

^{q) = e-'^ + e"'^ (29) 

This vector spans the one-dimensional physical Hilbert subspace Hp^yg C H . 

Let us check whether if) indeed projects functions onto their values at zero (up to 
the normalization constant). Let be any vector in H. Indeed, 

(^|</>) = f [(e-'^ + e'^-^) <i^{q) + (9, (e"" + e'^'^)) 9,</>(g)} (30) 

= \-e~'^ + e'^-') m]] (31) 

= c0(O) (32) 

with the normalization constant c = (1 — e~^). 

Finally, let us remark that since in this case Q is positive, and zero is therefore a 
boundary of the continuous spectrum, we could have found this solution also by using 
von Neumann's theory as mentioned above. 

Namely, the only boundary vectors, i.e. the only vectors in H obeying either 
(P* ± i)\ip) = 0, i.e. the only normalizable solutions obeying either of {—idq ± i)ip{q), 
are e^'^. Thus, we could have narrowed down the search for the physical subspace, 
knowing that it has to lie within the boundary space Hb spanned by these two vectors, 
as it does of course, being spanned by ip in Eq.p9[ 



5.2 Example: Similarly with isospinors 

Let us now consider an example of a constraint operator Q whose spectrum at zero is 
finitely degenerate, so that we can expect the corresponding physical, i.e. constrained 
Hilbert space to be multi-dimensional. 

To this end, let us consider the kinematics of a quantum mechanical particle which 
possesses isospin and which lives, say, on the positive half line. In its Hilbert space H 
the scalar product of wave functions then reads: 



N .00 
V^ / j^ J.* / ^\ J.I 



i=1^0 



') = E d^ ^K^)^'^{^) (33) 



We wish to constrain the particle from the bulk to its boundary by choosing a constraint 
operator Q, which acts as Qipi{x) = xipiix). By suitably imposing "Q = 0" we intend 
to project wave functions onto the A^-dimensional isospinor space at x = as the 
physical subspace. Within the original topology of H this is not possible because ker 
Q = ker Q* = 0. Following our general method, we therefore introduce the symmetric 
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operator P = —idx, obeying [Q, P] = il on its domain. We can now equip this domain, 
or more precisely, the domain Dp* , with the new scalar product Eq.p!7| 

AT 

to obtain the new Hilbert space H which possesses the improved topology. In H the 
operator Q is no longer self-adjoint and we identify Hphys = ker Q*. According to 
Eq.|l9|, the condition for vectors ?/) G -f/^ to be in the physical subspace, i.e. in ker Q*, 
now reads 

^ roo 

J2 / dx {0*(x)x^i(x) + id,x^*ix))dxiJiix)} =0 \f(l)eDQ (35) 

which yields 

x{dl - l)ipi{x) = (36) 

The solution space is spanned by the wave functions i/j^ (x) = 6n,ie~^, {n = 1, ...,iV), 
i.e. those vectors and their linear combinations represent Hp^ys in H. The projector 11 
onto the physical subspace is 11 = Yl,n=i IV'^"^) ® (V^*^"^]. 

Let us verify that the scalar product of the ^*^"^ with an arbitrary wave function 
(/)j(x) projects onto the isospinor space at zero. Indeed, 

(^W|0) = Y. rfx {(5„,,e-^-0,(x) + 9,e-XA0,(x)} (37) 

-e-^M^)]^ (38) 

= 0n(O) . (39) 

Clearly, our treatment of this example straightforwardly also applies in the case N = oo, 
to obtain an infinite dimensional physical sub Hilbert space spanned by the {c/)'^")}. This 
means that the classical constraint "manifold" which was here a set of A^ discrete points 
could also be taken to be an infinite set of discrete points. 

5.3 Example: Projecting onto the boundary of the half-plane 

Let us now consider the case where the classical constraint manifold is actually a 
continuous manifold. To this end we consider the example of a particle which lives, 
say, in a two-dimensional space on the half plane defined by Xi > 0, i.e. in its Hilbert 
space the scalar product of wave functions reads 

/oo TOO 

dx2 / dxi 0*(xi,X2)0'(xi,a;2). (40) 

-oo Jo 
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We choose Q = xi, in order to constrain the particle from the two-dimensional bulk 
to its one-dimensional boundary at xi = 0. Clearly, Q is positive, self-adjoint and 
possesses the half-axis xi > as its infinitely degenerate spectrum. 

To employ our method, we use the symmetric operator P := —idx^- As required, it 
obeys [(5,-P] = "il- Using P, the new scalar product reads ( | ) 



/OO /"OO 

dx2 / rfxi {r0'+(9,>*)9,,0'} (41) 

-OO JO 

which yields for our purposes the topologically improved Hilbert space H over Dp*. 
According to Eq.|19], physical vectors ij) G Hphys = ker Q* now obey 



/OO /'OO _ 

dx2 / dxi {^*xi^ + {d,,xi(l)*)d,,^} =0 y^EDQ, (42) 

-OO JO 

which means: 

Xi(92^-l)V^ = (43) 

Thus, the physical vectors ip ^ Hp^y^ are represented in H as wave functions ip{xi, X2) = 
f{x2)e~^^ where f{x2) is an arbitrary square integrable function. Let us explicitly 
verify that the scalar product of an arbitrary vector cf) E H with a physical vector 
ip G Hphys C H is the integral over the constraint surface. Indeed: 



00 


/•OO 


dX2 1 


-00 Jo 


00 p 


dx2 


— e 


-00 '- 


00 





(V^|0) = I__dx2 I dxi{f*{x2)e-'''^{xi,X2) + {d^J*{x2)e-''')d^,(p{xi,X2)} 
-e-^7*(a;2)0(xi,X2) 
dx2r{x2)(piO,X2) (44) 



00 




Q 

5.4 Example: Projecting onto a cylinder in IR 

For a less trivial example, let us now consider quantum mechanics in three space-time 
dimensions, choosing as the constraint operator: 

Q={X^,+Xl-R') (45) 

We then of course expect the constrained, physical Hilbert space to be the Hilbert 
space of square integrable functions over the cylinder with radius R around the x^ 
axis. 

In the original Hilbert space, H, the scalar product of wave functions reads, choosing 
cylindrical coordinates, 

/OO /'27r /"OO 
/ / dx3 dip dr r (PI cl)2 (46) 

-00 Jo Jo 
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while Q acts as 



Q 4>{r,^,X3) = (r^ 



R' 



[r,(p,X3) 



(47) 



As the first step, we introduce the operator P := —i^dr on its domain Dp of absolutely 
continuous square integrable functions (f) which vanish at r = and for which the almost 
everywhere defined derivative —i^dr4> is also square integrable. Second, we equip Dp* 
with the new scalar product, Eq.^, and its induced topology, to obtain the Hilbert 
space H. Third, we can now calculate the physical domain Hphys = ker Q*. According 
to Eq.|T^, a state ip is in the physical subspace exactly if for all e Dq it obeys: 

r27r roo 



/oo rzn /"OO / 

/ / dxs dip dr r (j)*{r'^ - R'^)ip + 

~oo JO Jo \ 

/oo /'2n /"OO / 

/ / dxs dip dr r (f)*{r^ - R^)%1) - 

-oo Jo Jo \ 



—dJr^-R^) 



drip 

2r ^. 



ir^ - R" 



^27r /"CxD 

+ / dip dxs 

J-oo 



^ir'-R') 



2r 



2r 

oo 




-dr 



2r 



dri) 



Thus, physical states must obey the boundary condition at r = 0, for all (p E Dq 

/•2tt /-oo ]^ ]^ 

lim / / dx^dip -(r^ -R^)(p*—drtp = 
r-->o+ Jo J-oo ^ ^ 2^ ^^ 2r ^ 



(48) 

(49) 
(50) 

(51) 



since the boundary term at infinity vanishes due to the square integrability. Also, for 
all positive r with r ^ R the physical states must obey the differential equation: 







r {r^ - R^)—dr—dr'ip 



(52) 



2r 2r 
The solutions to this differential equation are of the form: 

iP{r,ip,xs) = fiip,xs) {Ae'^-'+Be^') (53) 

Let us begin by considering the part of the solution in the region r < R: We know that 
ip must be continuous everywhere, and this is nontrivial at the origin, r = 0: Namely, 
for ip to be continuous at the origin one needs that either ip{r = 0, ip, x^) = 0, or that 
ip{r = 0,ip,Xs) is independent of ip. Thus, we obtain two different possible behaviors 
of physical states in the region r < R: The first set of physical states vanishes at the 
origin, i.e., 

V'(r,¥P,X3) = /(</^,X3) (e'"'-e"') . (54) 

Here, a priori, f{ip, x^) is some arbitrary square integrable function. However, solutions 
must also obey the boundary condition Eq.^ Since all (p either vanish at r = 0, or 
are yj-independent, Eq.^ takes the form: 



r2-K I 

= lim / dip — drip 

r^o+ Jo ^ 2r ^ 

r2iT 



lim / dip f{ip,X3) f-e "^ - e'' 

■^0+ Jo ^ 



(55) 
(56) 
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Thus, 

= -2 rdvf{^,xs) (57) 

Jo 

which means that in solutions of the form given in Eq.|5^ the /((/), X3) are indeed 
arbitrary square integrable function - with the exclusion of the y9-zero-modes. 

However, the zero-modes are not lost: The solutions which we have just obtained 
are only the set of solutions which are continuous at r = by virtue of vanishing there. 
As we said above, there is a second set of physical states, namely those which are 
continuous at r = 0. These do not depend on the variable ip, i.e. they are of the form 

ij{r,^,x,) = g{xs) {Ae-''" + Be'') . (58) 

For these solutions the boundary condition Eq.^ reads 

= lim f^dif —drip (59) 

= lim r d^g{x^) (-Ae-'^' + Be'") (60) 

r^0+ JO ^ ' 

= ~2tx g{x^){A-B) (61) 

which implies that B = A. Thus, the physical states which carry the zero-modes of 
the angular degree of freedom yj are represented in the region r < Rhy the functions: 

^(r,y,,a;3) = ^(0:3) (e-^' + e'"') (62) 

Clearly, for the part of the solution in the region r > R, square integrability requires 
solutions of the form of Eq.|5^ with B = 0. The solutions of the differential equation 
for r < R and for r > R are glued together by the requirement that every physical 
state is continuous at r = i?. We therefore obtain that the physical states are spanned 
by states of one of two forms. Either, 

(e''^ - e"""^), r < R 
ilj{r,(p,xs) = /((/?, X3) { (63) 



where / is an arbitrary square integrable function of X3 and ip, without, however, the 
zero-modes in (p, because of Eq.^. The rest of the physical states are the zero-modes 
in if, and are represented as functions of the form: 

[e-'-' + e''), r<R 
ip{r,Lp,X3) = gix3) { (64) 

3-"' (l + e^^'), r>R 

The physical Hilbert space Hp^ys C H is spanned by these functions, and is equipped 
with the scalar product of H, as given in Eq.p!7|. 
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In order to project an arbitrary function in H down to Hphys we need to take scalar 
products of arbitrary functions (p in H with functions in the physical Hilbert subspace 
Hphys- Let us check that, as desired, this scalar product reduces to an integral over the 
product of two functions over the surface of the cylinder r = R: 



(0,V^) 



/oo rZn /"OO / 

/ / dxs dip dr r ( 0*-?/' + 
oo JO Jo \ 

/oo r2n ( rR / ][ 

/ < / dxs d(f dr r ( (l)*ip — 0* — (9, 
~oo Jo \ Jo \ 2r 



2r 



drip 



+ 



+ 

2n POO 







oo /I 

dx3 dip dr r (jfip — (jf — dr 
R V 2r 

R- 



2r 

drip 

2r ^. 



dp) dx^ 



-d>* — drip 
2^ 2r ^. 



+ 



04 



2^ 2r 



drip 



R4 



(65) 
(66) 

(67) 
(68) 



where, to be precise, in line |6g the evaluations at the interval boundaries are limits 
taken from within the interval. The terms in lines ^ and ^ vanish because all physical 



states ip obey Eq.|5^. In order to evaluate line ^ we decompose (p into its ip- zero-modes 

00 and the rest, 0i: 

0(r,v9,X3) = (poir,X3) + (pi{r,p>,X3) (69) 

where 0i(r = 0,ip,X3) = 0. Recall that we found that while the physical subspace 
consists of all square integrable functions, say h{(p, X3) over X3 and over y^-space, its 
representation as functions in the original Hilbert space realizes the </)- zero-modes some- 
what specially. Namely, if we decompose an arbitrary h{p>, X3) into their (y9-zero-modes 
and the rest as 

h{^,X3) = f{ip,X3)+g{x3), 

then these are represented as 



(70) 



ip{r,(p,X3) 



We then read off: 



f{^,X3) {e-'-e~^')+g{x3Ke'-'+e-'''), r < R 



f{f,X3) e 



-r2 / 2ij2 



l)+gix3)e-''\l + e 



2i^2^ 



(71) 



r> R 



l^^rlpW 
l^dr1p\R_ 



2r 



dr1p\i 



V{f,X3) 

/(^,a;3)(e^'+e-«') + ^(0:3) (e^^ - e"^'') 
/(^,a;3)(-e^' + e-^') - ^(a;3)(e^' + e"^') 



Using 



0i(O,v2,a;3) 



0, 



(72) 
(73) 
(74) 

(75) 
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dip (l)l{r,(p,xs)g{xs) = (76) 

/o 

and 

r"rf^05(r,a;3)/(<^,X3) = (77) 

Jo 

we therefore find that: 

(0)^) = 27 y^^^V'c^Xs |(0o(r,X3) +0i(r,(/?,X3))— (9^^|,.=fl_ (78) 

- i<P*oir, X3) + 0*(r, If, X3)) —dr'il>\r=o+ (79) 

- (0S(r,X3)+0t(r,^,X3))^9,7/>|,=H+} (80) 

„ ("277 ("OO 

= e^ / / rf<^dx3 {0S(i?,X3)^(x3) + 0t(^,<^,^3)/(<^,a;3)} (81) 

2 /'27r /"OO 

= e^ / / rf¥Pc?X3 0*(/?,¥^,X3)/i(<^,X3) (82) 

JO J-00 

as it should be, with the factor e^ being an overall normalization constant. 

5.5 Example: Projecting onto a line in IR 

Let us briefly also consider the case of projecting not onto a cylinder, but onto a 
one-dimensional line in IR^. To this end, we consider, similar to the previous example: 

Q = xj + xl (83) 

and again P = —i^dr, where r^ = x^ + x^. Within our formalism, the change of 
the dimensionality of the constraint manifold is automatically taken care of: The dif- 
ferential equation obeyed by the physical states, Eq.^ must now hold for all pos- 
itive r > 0. Thus, normalizability shows that all solutions must be of the form: 
il!{r,ip,X3) = f{(f,X3)e~^ . Since these functions do not vanish at the origin, r = 0, 
continuity at r = requires independence of ip: 

ip{r,ip,X3) = g{x3)e~'' (84) 

While for R> these solutions were ruled out by the boundary condition Eq.^, here 
Eq.^ is trivially obeyed due to its prefactor r. 

Let us check whether the scalar product of any physical state ip with an arbitrary 
state reduces to the integral of the product of the two functions over the X3 axis. 
Indeed, 



00 /•27r 



(01^) = / dx3 dip 
-00 Jo 

(•2-K \ 



-(p*{r,^,x-i)—dri) 



(85) 

/OO /"ZTT I 

rfxg / dip -(j)*{Q,ip,X3)g{x3) (86) 

-00 Jo I 

/OO 
rfX3 0*(O,X3)^(x3) (87) 

-00 
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where in the last step we used that the functions (f) do not depend on y? at r = 0. 
Thus, the scalar product of arbitrary functions in H with a function from the physical 
subspace indeed yields the integral of the product of the two functions over the X3 axis, 
as expected. 

This example also shows how the formalism takes care of a dimensional reduction: 
The loss of the (/^-degree of freedom is expressed by the representatives of physical states 
possessing only zero-modes of the (y9-degree of freedom. In other words, restricting the 
function space over a cylinder to only its zero-modes in the angle (p is the same as 
reducing the function space to a function space over a line rather than the cylinder. 

5.6 Example: An impossible constraint 

It is perhaps instructive to briefly discuss how an "impossible" constraint yields an 
empty physical subspace. To this end, let us consider again in IR^ the constraint 

Q = xj + xl + R^ (88) 

for some positive R > 0. In this case, physical states must obey the differential equation 



Eq.^ for all positive r > and therefore, by the above arguments, the only possibility 
are functions of the form given in Eq.Q. However, unlike in the case R = 0, now the 
boundary condition Eq.^is nontrivial, namely yielding g{x3) = 0. This rules out any 
solutions, thereby yielding Hphys = 0, as should be the case. 

6 Uniqueness of Q and P, and the topology of the 
constraint manifold 

As we saw already in Sec.|^, our method cannot be unique: While the topology of the 
Sobolev space H^ is good enough to allow the projection onto a function's value at 
a point, higher Sobolev spaces H"- could be used as well. Analogously, in our more 
general method, for a given choice of Q and P, the topology of H allows us to project 
onto the physical Hilbert subspace, while higher Sobolev space analogs of H could be 
used as well. Interestingly, however, there is even more non-uniqueness in our general 
method. 

Firstly, for a given classical constraint manifold there is a non-uniqueness in the 
possible choices of the constraint operator Q. Secondly, also for a given Q there exists 
in general a non-uniqueness in the possible choices of a canonical conjugate operator 
P. The question in how far the various physical Hilbert subspaces obtained from these 
different choices are or are not equivalent is subtle and reveals interesting relations to 
the dimensionality and topology of the constraint manifold. 
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6.1 Nonuniqueness of Q for a given classical constraint 

Perhaps surprisingly, even simple classically identical constraints may differ quantum 
mechanically, even when there are no ordering ambiguities: In order to see this, let us 
consider the example of the constraint operator Qi = x\ + x\ + x\ — R^ in the Hilbert 
space H of square integrable functions in IR^. Classically, the constraint manifold, i.e. 
the set of solutions to Qi = is a 2-sphere in IR^: The physical Hilbert space which we 
obtain from Qi is of course the L^ Hilbert space of functions over this sphere. Similarly, 
the constraint operator Q2 = x\ + {x\ + x^ — R^Y describes the equator of the sphere 
described by Qi. 

Now, let us consider the constraint operator Q^ := QiQ2- Classically, its constraint 
manifold, i.e. the set of points for which Qs = 0, is still the 2-sphere with radius R, i.e. 
it is identical to the case above. Quantum mechanically, the physical Hilbert space for 
Q2, however, can now be different from the space of functions over the sphere. This 
is because Q2 counts the 1-sphere of the equator twice. We might expect that the 
physical Hilbert space also specifically contains L^ functions over the equator and their 
scalar product. In this way, the scalar product could be consisting of evaluating both 
the two-dimensional integral over the sphere and the one- dimensional integral over the 
equator of the product of two wave functions. 

This issue will be further investigated elsewhere. It is clear that whenever such 
ambiguities arise, i.e. whenever the physical subspaces can be, for example, function 
spaces over a variety of topologically different manifolds (of varying, and in general 
possibly not even well-defined dimensionality), then this should also have a reflection 
in the functional analysis of P. 

6.2 Non- uniqueness of the choice of P for a given Q 

For a given choice of Q, our condition that P be a maximal symmetric operator obeying 
[Q, P] = il does not determine the choice of operator P uniquely. Clearly, any two 
operators P and P' which obey the commutation relation [Q, P] = il will differ by a 
symmetric operator R which commutes with Q, i.e. P' = P + R with [R, Q] = 0, with 
the appropriate domains understood. Usually, two such operators P will be connected 
by a gauge transformation, in the sense that there exists a unitary operator U in H, 
such that P' = P + R = UPW . In this case, as is easily verified, U is also an isometry 
connecting the corresponding "generalized Sobolev spaces" H and H', as well as being 
an isometry connecting the corresponding physical subspaces Hphys and H'^^y^- In this 
case, the choice of P is therefore immaterial. 

For example, if two operators P differ by an operator R which is a self-adjoint 
operator that is independent of the constraint degree of freedom, in the sense that it 
commutes with Q and P, then we find U = exp{iQR). Or, if P is a function k{Q) of 
Q with integral J dq k{q) = K{q), then we have U = exp{iK{Q)). 

More concretely, let us reconsider the example of Sec.^.3|, where we used Q = Xi to 
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project onto the boundary (the X2-axis) of the half-plane {xi > 0). We had chosen P = 
—idx^- But this choice was not unique because obviously e.g. also P' := —idx^ + g{x2) 
obeys [Q^P'] = il. The question arises as to whether we obtain the same physical 
Hilbert space if we use P' rather than P. A short calculation shows that, using P', 
the physical vectors are represented as functions f{x2)e~^'^^^~^^^^^^^\ The difference 
therefore is only a phase, which vanishes for xi = 0, i.e. after the projection. Thus, 
the use of P' yields the same Hilbert space as does the use of P. 

Similarly, also e.g. P" := —idx^ — idx2 obeys the commutation relation [Q, P] = «1. 
Using P", physical vectors must obey the differential equation {dx^ + dx^ — I)'?/' = 0. 
Using the coordinates z := {xi + x-i)l2 and y := {x\ — X2)/2, we obtain (pz — 1)ip = 0, 
yielding for the general expression of physical vectors: ip = f{—y)e~~^, i.e. ip = f{x2 — 
Xi)e~^^^~^^'^\ The scalar product of a general vector (p E H with a physical vector 
^ e Hphys is then: 

rl^. I ^^. fj,*^J. {(f) J. f) \f*f-^.—-^J\^-i^l+^2)/A(fi J- f 

X2 





OO 

dx', 

-OO 
OO 



-OO 
OO 



{^\ct>) = 1^ dx, l_^dx2{r<P+{{dx,+dx,)nx,-X2)e-'^'-''+'-'''>/'){dx,+dx 

dX2^p*iO,X2)(l)iO,X2) (89) 



Thus, when using P" we still obtain the same physical Hilbert subspace. 

Similarly, also P'" := —idx^ + g{xi) obeys the commutation relation. Its use leads 
to physical vectors being represented as functions ip = f{x2)e~^^~^J amd^:^ ^jj ^g 
obtain is a new phase factor which vanishes (up to possibly an irrelevant global overall 
phase) after the projection onto Xi = 0. Thus, we obtain the same Hilbert space also 
from P'", as long as g{xi) is integrable. 

Even though these examples illustrate the robustness of the projection onto the 
physical subspace under changes of P we cannot exclude that, in general, R may 
not be a pure gauge, i.e. it may not be "integrable" in these simple ways. One 
may find unitarily non-equivalent maximal representations of the commutation relation 
[Q, -P] = il. On the other hand, while a priori these representations will lead to non- 
equivalent generalized Sobolev spaces, some of the resulting representations on the 
corresponding physical subspaces may still be equivalent after projection. 

The investigation of the precise relation between different choices of Q for the 
same classical constraint, unitarily non-equivalent P's, and the non-uniqueness of the 
topology of the "constraint manifold" is probably deep and should be worth pursuing 
further. 
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